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INTRODUCTION 


This report details results of research performed during the period April 1, 1989 through 
March 31, 1990 under NASA contract no. NAG-1-979. 

The general objective of this research has been to develop a life prediction methodology 
for laminated continuous fiber composites subjected to fatigue loading conditions. 

SUMMARY OF COMPLETED RESEARCH 

The following is a summary of research completed during the preceeding contract year: 

1) a phenomenological damage evolution law has been formulated for matrix cracking 
which is independent of stacking sequence; 

2) mechanistic and physical support has been developed for the phenomenological 
evolution law proposed in 1); 

3) the damage evolution law proposed in 1) has been implemented to a finite element 
computer program; and 

4) preliminary predictions have been obtained for a structural component undergoing 
fatigue loading induced damage. 

The results reported above are covered in detail in the enclosed appendices. 
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Abstract 

The accumulation of matrix cracking is examined using continuum damage 
mechanics lamination theory. A phenomenologically based damage e\o 
lutionarv relationship is proposed for matrix cracking in continuous fiber 
reinforced laminated composites. The use of material dependent properties 
and damage dependent laminate averaged ply stresses in this evolutionary 
relationship permits its application independently of the laminate stacking 
sequence. Several load histories are applied tocrossply laminates using this 
model and the results are compared to published experimental data. The 
stress redistribution among the plies during the accumulation of matrix 
damage is also examined. It is concluded that characteristics of the stress 
redistribution process could assist in the analysis of the progressive failure 
process in laminated composites. 


Introduction 

A unique propertv of composite materials is their evolutionary failure char- 
acteristic. The inhomogeneity of the microstructure provides numerous 
paths in which loads can be redistributed around the damaged region. 
Thus, the integrity and response of the component are affected more by 
the collective effects of the accumulated subcritical damage than by any 
single damage event. For laminated composites, this subcritical damage 
takes the form of matrix cracks, delaminations, debonding, and fiber frac- 
tures. Because the transfer of load away from the damaged area influences 
the damage evolution in the adjacent areas, the stresses at the ply level play 
an important role in the evolution of damage and the ultimate failure of 
the structure. Thus, this unique failure process that makes composites an 
attractive engineering material has also limited their efficient use in struc- 
tures. The primary reason has been the shortage of analytical means to 




mode] the complex events occurring within the laminate and the prediction 
of damage evolution. 

A review of the current literature indicates that three approaches have 
been taken to solve this problem. The first approach is the phenomeno- 
iogica 5 methodology. Empirical relationship" are developed from a body o: 
experimentally measured data. Statistical theory is frequently employed to 
correlate this data. Models of this type teno to be restricted to the stackim 
sequence utilized to construct them. The second approach., called the crac*. 
propagation method, identifies damage as dominant cracks and fracture 
mechanics is applied to predict crack growth. The physical significance of 
each damage mode is retained with this approach. Unfortunately, the dam- 
age state in composite materials contains a multitude of interacting defects. 
Analysis in this manner is complex and perhaps unmanageable. For exam- 
ple. if there are hundreds of cracks, then the hniie element method .would 
require lens of thousands of elements. The third and most recent approach 
is the use of internal state variables in a continuum damage mechanics 
framework to model the damage accumulation. The averaged effects of the 
damage are represented through the internal slate variables. This theory 
provides a thermodynamically rigorous characterization of the continuum 
under examination. The continuum damage mechanics approach entails 
the formulation of constitutive relationships, damage variable descriptions, 
and the damage evolution laws. The damage evolution laws may either be 
phenomenological, mechanistic, or even some combination of both. Because 
this approach is capable of accounting for the slacking sequence, it differs 
considerably from the phenomenological approach. Furthermore, because 
the effects of each crack are treated in the constitutive equations rather 
than via fracture mechanics, t he computational solution is simpler than 
the crack propagation approach. 

In the current paper, the continuum damage mechanics lamination the- 
ory proposed by Allen, e; al. (l9S7a.b) is reviewed. Also a matrix crack 
damage evolutionary law is developed and this phenomenologically based 
evolutionary relationship is used to examine the accumulation of damage 
and the accompanying stress redistribution among the plies ir. laminates 
subjected to fatigue loading conditions. 


Review of the Damage Dependent Lamination Model 

The continuum damage mechanics approach used herein is based on the 
thermodynamics of irreversible processes. It is postulated that the state of 
a materia] point in a system undergoing a dissipative process can be char* 
anerized by a set of observable and interna] state variables if this process is 
sufficiently close to the equilibrium state. Through the application of con- 
straints from the fundamental principles of thermodynamics as proposed by 
Coleman and Mize! (1966) and the assumption of interdependence among 
these state variables. Coleman and Gurtir. (1967) showed that constitutive 




equations for the material can be constructed in terms of the strain tensor, 
temperature, and internal state variables. The internal state variables may 
represent any dissipative process occurring in the medium. For the current 
application to distributed matrix cracks ir, laminated polymeric compos- 
ites. Alicr. e*. a!. (3 9S7a'j selected a second order tens'?: value interna: state 
variable to represent the kinematics of the cracks. This tensor was first 
defined by Vakulenko and Kachanov (1971) to be 


Qi; = — [ u x n : dS, 
1 L JS 
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where a., are the components of the internal state variable tensor. 1 i is 
a local volume in which statistical homogeneity can be assumed, it. and 
7 i. are the crack face displacement and normals, respectively, and 5 is the 
crack surface area. 

The thermomechanical response of an elastic material with damage was 
found to be related to the damage dependent Helmholtz free energy as 
follows, 


G Li j 


dhi 
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where are the components of the locally averaged stress tensor, hi is 
the volume averaged Helmholtz free energy and r i, : are components of the 
infinitesimal locally averaged strain tensor. The subscript L will be used 
herein to denote volume averaged quantities. The Helmholtz free energy, 
hi, for an elastic material with distributed damage is denned by 


hi = hv L - u e L% 


( 3 ) 


where hri is the Helmholtz free energy of an equivalent undamaged elastic 
body ana u c L is the energy associated with the damage effect on the ma- 
terial. These two energy quantities are expressed by second order Taylor 
series expansions of the corresponding state variables. If the higher ordered 
and residual effects are neglected and isothermal conditions are assumed, 
The ply level thermomechanical constitutive relationship of the damaged 
material is, 

C\j — CnjklZkl -r iLijki&kl- ( 4 ) 

where Ci-jki is the undamaged material modulus tensor and li,jki * s de- 
noted the damage modulus tensor. The components of this damage modu- 
lus tensor were shown by Allen, et al. (1987b) to be related to the modulus 
tensor as follows 

?Lijkl ^ —Clijkl- (^) 

Another result of the Coleman-Mizel formulation is the entropy production 
by the distributed damage in the absence of thermal gradients. 




where the thermodynamic force, jj } , associated with the evolution of o i; is 
denoted by 


inequality (6) will admit only those processes that yield non-negative rates 
of entropy production. 

The response of a multilayered laminate with matrix damage is obtained 
by inplane averaging the ply level constitutive relationship show in equa- 
tion (4) and imposing the KirchhofT hypothesis through the thickness to 
obtain the following modified set of laminate equations. 
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where {.V} is the resultant force per unit length vector, [Q jfc is the elastic 
modulus matrix for the J: th ply in laminate coordinates, a* is the distance 
from the midplane to the k ihi ply, c° and k are the midplane strains and 
curvatures, respectively; and } are the damage variables corresponding 
to the matrix cracks in the J: tr ' ply. The effects of the matrix cracking on the 
resultant forces are contained in the last group of terms on the right hand 
side of equation (6). An expression for the resultant moments can be written 
in a similar form. The stiffness loss of composite crossply laminates with 
matrix cracks has been successfully predicted with this model by Allen, 
et a!. (19S7b, 1986) for specified damage states using equation (8). In 

addition, although it will not be discussed in this paper, the model has also 
been extended to include the effects of delamination by Harris, et al. (IPS* , 
198S). A damage evolution relationship for this model will be introduced 
in the following section. 


Evolutionary Equation for Internal State Variables 

In order to characterize the development of the internal damage, the con- 
ditions for the initiation and progression of damage as well as a means to 
quantify changes in the damage state must be established. Concepts from 
linear elastic fracture mechanics can often be employed to assist in the 
development of the evolutionary equations for brittle damage. The ther- 
modynamic force, f?. as defined by equation (7), represents the available 
free energy that can be delivered by the system for a small change in the 
internal state variable. However, this change will occur only if the energy 
delivered is equal to or greater than the energy required to produce this 
change in the internal slate. Thus, the following initiation criterion based 




on the relative magnitudes of the available and required thermodynamic 
forces. (ijA can be used for this model, 

' 1 J 1 r e q. 



The implementation of this criterion will require the determination of the 
required thermodynamic forces. The required thermodynamic forces are 
most likely to change with the damage state. This criterion is analogous 
to comparing the strain energy release rate to the critical value in fracture 
mechanics. Had the interna! state variable been defined as the crack sur- 
face area, the thermodynamic force would be identical to the strain energy 
release rale. 

One approach to the formulation of the internal state variable evolution- 
ary relationships is through micromechanical considerations. However, this 
approach is depend cm on the availability of micromechanical solutions that 
can model the essential physical characteristics of the damage state. For 
the problem of matrix cracks embedded in an orthotropic medium that is 
layered between two other orthotropic media, the solutions that are cur- 
rently available are applicable only to very specific loading conditions and 
damage geometries. Therefore, the evolutionary equation proposed herein 
is phenomenological in nature. The form of the damage evolutionary rela- 
tionship employed in this paper is based on the observation made by Wang, 
ei al. (198-4) that lor some materials the rate of damage surface evolution 
per load step. follows the power law as shown below, in which the 

strain energy release rate, G, and a material parameter, n, serves as the 
basis and exponent, respectively. 


^ = PG r - (10) 

where P is a material constant. To develop internal state variable evolu- 
tionary equations in the form of equation (10). a :j has to be related to 
the damage surface area. Since a j; , as defined by equation (1), represents 
the kinematics of the crack faces, the damage surface area alone will not 
be sufficient to describe the crack face displacements. Assuming that each 
crack in the material volume shares a common geometry and orientation, 
then the specification of the far field strains will complete this description. 
The rate of change of the internal state variable can therefore be expressed 
as follows. 

_ QQ U ^ 

dN ~ dS dN 1 ; 

where the far field strains are reflected in the term which relates the 

d .> * 

changes in the internal state variable during the development of damage 
surfaces. Thus, using equations (10) and ( 1 1 ) , the stable evolution of the 
internal state variable is defined by 
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where can be determined analytically for simplified damage geometries 
and loading conditions and k j is a material parameter. 

Jl will now be shown that the strain energy release rate can be determined 
from the thermodynamic forces. The Helmholtz free energv for a body with 
distributed damage was defined earlier to be the sum of the Helmholtz free 
energv in an equivalent undamage body, and the energ> of damage 

creation, Since h£i is independent of the internal state variables, the 
thermodynamic force can be expressed in terms of as 





(13) 


Allen, el al. (1987a) defined u' L as the mechanical energy of a continuum 
due to equivalent crack surface tractions acting on the crack faces, i nis 
energy encompasses the energy available for crack extension an C me energ} 
loss due to the apparent stiffness reduction caused by existing cracKS. --m 
association between the energy for crack extension and the strain energy 
release rate, G. is defined by 


1 

vL 



(14) 


The strain energv release rate is obtained from equation (14) by different! 
ation with respect to the crack surface area. Thus, 


cu) 

G = V L — 
d$ 


(15) 


If the process is restricted to isothermal conditions, v c L depends only on the 
strain tensor, c, J? and the internal state variable. a tJ , so that equation (15) 
becomes 

(16) 
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The relationship between the strain energy release rate and the available 
thermodynamic force is obtained by permitting crack extension under fixed 
grip conditions and using the definition of the thermodynamic force in 
equation (13). 

r _ yd*»,T ( 17 ) 

G ~ L-rz-Jxy ^ ' 

The internal state variable evolutionary equation shown in (12) can thus 
be expressed in terms of the available thermodynamic force as follows: 
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when f* > (/,;),„ • ( 18 ) 


For crack propagation under a single fracture mode, the crack surface kine- 
matics for a thin laminate can be characterized by a single component of 




the internal state variable tensor. For example, consider the matrix crack 
damage state shown in Fig. 1, where all the crack planes are fiat and ori- 
ented perpendicular to the plane formed by the lamina and parallel to the 
fibers, a pure opening mode (mode I) would correspond io the crack faces 
moving in a direction paraiiei to tne crack face normals. Tnus, in reference 
to the ply level cartesian coordinate system, the internal state variable ten- 
sor as defined by equation (1) will be zero except for 022 - Likewise, for the 
pure shearing mode (mode II), the only non-zero component of the internal 
state variable tensor would be o;j. 


Matrix Crack Development in Multi-Ply Laminates 

The evolutionary relaiionsnip shown in equation (16) provides a description 
01 the damage development at the ply level. In laminates containing multi- 
ple plies of different orientation, the development of damage is influenced by 
the adjacent plies. This is attributed to the different stress states found in 
the neighboring plies and the redistribution of load among the plies that oc- 
curs during the accumulation of damage. The effects of this interaction on 
the damage development are included implicitly in the evolutionary equa- 
tion through the thermodynamic forces and because laminate averaged 
damage dependent ply responses calculated from equation (8) are used to 
determine these quantities. Since the adjacent ply constraining effects are 
reflected through the laminate equations, the evoluti^Ul' equation can be 
applied independent of the laminate stacking sequence. 



Figure 1. localized configuration of ply with matrix cracking. 
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The term , found in equation (18), reflects the changes to the inter- 
nal state variable with respect to changes to the damage surfaces. can 

be obtained analytically from relationships describing the kinematics of the 
crack surfaces for given damage states and loading conditions, should such 
solutions exist. For transverse matrix cracks in crosspiy laminates, tne aver- 
age crack face displacement in the pure opening mode can be approximated 
by a solution obtained by Lee, et al. (1989) for a medium containing an in- 
finite number of alternating 0 C and 90° plies. Thus can be determined 
for crosspiy laminates subjected to uniaxial loading conditions. It has been 
found that for typical continuous fiber reinforced Graphite/Epoxy systems 
can be assumed to be constant for a given applied load until the dam- 
age state has reached an advanced stage of development. This assumption 
has facilitated the determination of the material parameters /;j anc n. The 
damage state at any point in the loading history can now be determined by 
the integration of equation (IS) using the laminate averaged ply responses. 
This integration is performed numerically because of the nonlinearity of the 
damage evolutionary equation. The fourth order R.unge-Kutta method has 
been found to be suitable for this application. 

The development of the matrix crack damage state in crosspiy laminates 
subjected to uniaxial fatigue loading is examined using the proposed dam- 
age evolution equation. To maintain the thermodynamic admissibility of 
the fatigue loading process, it is assumed that the values of the internal 
state variables remain constant during the unloading portion of the load 
cycle. It was further assumed that the required thermodynamic force is very 
small compared to the available thermodynamic force, thus a 22 will change 
at the onset of load application. The material properties for AS4/3503- 
6 Graphite/Epoxy are used in the calculations to enable the comparison 
of model prediction 10 experimental measurements made by Chou, et al. 
(1952). The material parameters for this polymeric composite system have 
been found to be 

k\ = 4.42 and n = 6.39 . (19) 

The damage histories for two crosspiy layups have been predicted using the 
model. One layup consists of four consecutive 90° plies laminated between 
0° layers of two ply thickness. The other crosspiy laminate contains six con- 
secutive 90° plies in the center. The model predictions of the damage state 
in the 1 O 2 / 9 O 2 I 5 laminates fatigue loaded at a maximum stress amplitude of 
38 ksi and 43 ksi are shown in Figs. 2 and 3, respectively. The lower stress 
amplitude is equivalent to eighty percent of the monotonic crack initiation 
stress, while the higher stress amplitude is equal to ninety percent of the ini- 
tiation stress. The experimentally measured damage states were originally 
measured in terms of the crack density. However, the corresponding 022 
for each damage level can be approximated by the relationship proposed 
by Lee, et al. (19S9). 









Good agreement is found between the model predictions and the experi- 
mental results. The damage evolution for the thicker |0 2 /903], laminate is 
shown in Fig. 4. This laminate was loaded at a maximum stress amplitude 
of 26 ksi. This amplitude corresponds to eighty percent of the quasi-static 
matrix crack initiation stress. The results for this load case indicated pood 
agreement with the experimental data. The effect of the load redistribution 
or the damage evolution is apparent in this load case. A marked decrease 
in the rate of damage evolution after fifty thousand load cycles was indi- 
cated by the model. On the other hand, the experimental data showed this 
decrease to occur after only ten thousand load cycles. Since the evolution of 
delaminations were not included in the analysis, the predicted decrease in 
the rate of damage evolution is attributed only to the matrix crack induced 
transfer of load from the 90° plies to the adjacent 0° plies and the resulting 
decrease in the available thermodynamic force. The measured values of 
the damage state, however, may have been influenced by the formation of 
delaminations along the free edges and in the interior. Such occurrence can 
ciastically affect the stress distribution among the plies and the available 
thermodynamic force for damage evolution. 

To examine the amount of stress redistribution that occurs during the 
damage accumulation, the model was used to determine the axial stress in 
the 90° plies of the [ 02 / 903 ], laminate fatigue loaded at three different stress 
amplitudes. Fig. 5 shows that for the stress amplitude of 38 ksi. the axial 
stress in the 90° plies after forty thousand cycles was less than fitly per- 
cent of the original stress level in the undamaged laminate. Therefore, the 
rate of damage evolution is expected to be relatively low during the latter 
stages of the loading history. This is observed in Fig. 6. which shows the 
corresponding values of the internal state variable, 022 ? for the three loaa 
cases. The 26 ksi stress amplitude load case, on the other hand, produced 
only gradual changes in the axial stress and damage state as compared to 
the other two stress amplitudes. The percentage decrease from the original 
undamage stress level increased with the fatigue stress amplitude. These 
results demonstrate that the stress redistribution characteristics among the 
plies in the laminate are dependent on the loading conditions. These redis- 
tribution characteristics will affect the manner in which damage develops 
in the surrounding plies as well as eventual failure of the laminate. 


Conclusion 

A damage evolutionary relationship for the accumulation of matrix cracks 
has been presented. This phenomenologically based relationship operates 
as part of a continuum damage mechanics model developed to analyze the 
response of laminated composites. The utilization of material dependent 
quantities and damage dependent laminate averaged ply responses in the 
evolutionary relationship has enable it to function independent of the lam- 





Figure 4. Matrix crack damage in the 90 p plies of a [0^/9 O 3 j , AS4/3501-6 
laminate loaded at a stress amplitude of 2 oksi . 



Figure 5. Damage induced stress redistribution in the 90 c plies of [O 2 / 9 O 3 ], 
laminates subjected to constant stress amplitude fatigue loading. 







Figure 6 . The accumulation of damage in the 90° plies of IO 2 . y 9 0 3 u «lamin 
-ales subjected 10 constant stress amplitude fatigue loading. 

inate stacking sequence. The capability to predict the evolution of matrix 
cracks in crossply laminates subjected to fatigue loading conditions has 
been demonstrated. The evolutionary relationship has also been used to 
examined stress redistribution among the plies during the damage history of 
the laminate. It is found that the stress distribution behavior is dependent 
on the load amplitude. Higher applied loads result in rapid changes in 
the axial stress and damage state during the initial portion of the loading 
history. This is then followed by low rates of change for the remainder of 
the loading history. For lower load amplitudes, the stress redistribution 
process and damage accumulation proceed in a gradual manner. 

The existence of other types of damage, such as delaminations, will alter 
the redistribution of stress among the plies and thus the evolution of the 
damage state. An evolutionary relationship for delamination is currently 
under development. The inclusion of this type of damage should provide 
better insight into the complex events occurring within the laminate. The 
information obtained from the stress redistribution histories could enhance 
the knowledge of failure characteristics in laminated composites and aid in 
the development of models for this progressive failure process. 
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A MODEL FOR PREDICTING DAMAGE INDUCED fATIGUE LIFE 
OF LAMINATED COMPOSITE STRUCTURAL COMPONENTS 
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Tr.ii paper present- a model for predicting 
tne life or laminated composite structural 
components subjected to tatigue induced 
mi crost rue tuna 1 carriage. Tne mooe ' uses tne 
concept of continuum odmace mecnanics, wherein 
the effects of microcrac». s are incorporated into 
a damage dependent lamination theory instead of 
treat i no each crack as an internal boundary. 
Internal variables are formulated to account for 
the effects of both matri> cracks and internal 
de 1 ami nat i ons . Evolution lews for determining 
tne damage variables as functions of ply stresses 
are proposed, and corr.pari sons of predicted damage 
evolution are made to experiment. In addition, 
predicted stiffness losses, as well as ply 
stresses are snown as functions or damage state 
for a variety of stacking sequences. 
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Structural applications using composi^r 
materials continue to increase in tne aerospace 
engineering community. 1 ms *S primarily due to 
me fact inat aerospace applications are driven 
by Mer. suecmic strencm accompanied by low 
mass. Altnough compos ilfrS are expensive wnen 
compared to monolithic materials, their improved 
properties can make them cost effective 
alternatives in suer, applications as soace 
structures, therein tne cost of launching a pound 
of mass is exorbitant. 

A*':noucn composites possess many inviting 
material cr.aracter i sties, mere are nevertheless 
some sno*‘tcomings. Possibly tne most sic^incant 
one is tne : ac : tnai most composites ceveiop 
mi crosiruciura 1 carnage when subjected to long 
term fatigue condition-. Tnis is cue to tne 'act 
that combining two materials with dissimilar 
mecnanical properties results in a multitude or 
stress singularities when external ioacs are 
applied. These tremendous gradients contribute 
to stable nncrocraking ».r iCh continues to evolve 
tnrougnout Hit life of ire structural component, 
and, in mcM ci rcumstanceu , tne damage leaas to 
ultimate fracture of tne part ana subsequent 
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Figs. 1-4 i.*), wherein successive damage stales 
are shown ror a coupon of graphite/epoxy (0/90^ ) s 
composite Subjected to ratigue loading. In the 
photos the coupons have been subjected to 
uniaxial loading in the vertical direction, 
leading to transverse matrix cracks in the 9CT 
plies (Horizontal lines), longitudinal splits in 
the 0 piles (vertical lines) and interply 

de lamina t ions between the IT and 90' plies 
( s ho Jed gru, areas). 1 p 1 Cal 1> , the n>zlr U 
craCH* induce stress concent rat ions wnich preceor 

t f »~ e : 0 . L' jr i an* i f.c U*' k . at Sh:^r ; ir r , j. . 

5. i u ri heriitorc , inis jaiuntjc is accompanied Uj c 
small toss m component c > i c stiffness, as snown 
ir, Fig. i. it :,cs oeer. conclusively determined 
that tne carnage is stress induced, so tnat the 
damage state and resulting stiffness are strongly 
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add it ion. 

the 

growth 

CT 

damage is 

dependent 

on stac 

k i ng 

sequence 

, c : 

i 1 1 ustrat cd 

i r, Fic. 

? , hot e 

tnat 

the dele 

in ha 

cr.s tend to 

propagate 

normal tc 

tne 

lea: in: 

C i 

c : t ion in a 

(0/90 j) s 

laminate , 

as 

cp pose: 

to . 

k*r; ical ly in <■ 

: (0/90, )< 

laminate. 


is r.ypot 

ne i 

::cc that tnes 

e variitii 

3ns are du 

e :: 

differences 

i r. ply 1 e . c 

• ! stresses caused 

by 

chanoi no 

1 1 

p 1 j St 6Ck 1 n 

g sequent 

:e . In f 

act , 


this process can be considered to be in some 
sense beneficial because the integration of two 
separate brittle materials results in damage, 
which is o ductile-like phenomenon, thus causing 
load transrerral not unlike that which occurs in 
metals. 

The current procedure for predicting life m 
compos lies seems to be largely phenomenological 
in nature; that is, the metnua of analysis IS in 
most casei dependent or .opponent ceome:^. lea: 
histcr,, ever, vexing sequence an: 

fn v 1 r O f VT , e'* t z C Or.*. 1 1 1 1 C ’ ■_ ^rCDdt “ t rc rTl*. 

direct ;r c w l cure utilizer ad noc failure lc*i net 
unhxt Miner* c r_ie ! j , j ; . As a result, tr, C ; 
apply tin!, a preoc i c rr. i ned geometry an: 

stacking ieqjence. F'cri.ai.c tne most ambitious 
aoprcacr. use^ linear elastic finite elements an: 
treats eacn era ck as an internal boundary subject 
tc growth ueterminec fracture mecnanics 

criieric | 4 , 5 i. Wni k- this is desirable and 
pr obdb 1 j t i.e most accurate approach, it «Oulc 
require supercomput ing capability for a typically 
complicated damage state. 

For the pail Severe* years the authors hcv- 
been developing an alternative model for 
predicting damage development in laminated 
compos ik : . i ne moo el can be used for any 



Fie. 1, [marge X-ray radiograph or r 1 { O/ 9 G 2 > s Fie. 2. tnla^oed X-ray radiograph c: *1 (G/yO^/c 

f =. 2 . 0 n z , k = C . 1 , 


.v£, 0 = 250,00:, {' = 2.0 nz, k=Cj.i, 

. n _ i 10.230 km 


= 1C. 233 xt. 
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stacking sequence, using only initial linear 
elastic ortnctropic ply properties. It has the 
capability u. predict t tie effects of matrix 

cracks <ino oe lami nat i ons on ply stresses. 
Therefore, n can De useo to predict damage 

evolution c: a function of load history such as 
that snown i.\ } igs. \-A. I ur thermo re , the mode i 
fliaj be UiL-t -w mouc'i 1 1« cr responses Of bear s 0 nj 

plates v.ur stress gradients, so that it is not 

component . • m c. /■: me time of inis witirc 

trie moje • i.a. reached „ t c i ^ 1 > acver.ced State, 
a 1 tnc-q* ' u r>r ; ■ r .i.j : .me * c c ; 1 ' me mode 
nature ^ i-e'nt. ''■!£ Piper give: a o- u * 

r e » i r » • : i'*c lutsp: si etc cr oeve s cement of ;ne 

mccei. ’ ’ r — , c; r ».e j. nc: pretend to Suggest 
that t n- Kiooe ! is a usable cesign tod a: this 
time, v.e ere hopeful that the methodology 
suggested nerein is iiu.ing in the right 

di rec: ion. 


Hope! Oeve ippment 

T ne structure is muoelled as a simply 
connected domain, *itn the effects of 
microcrac* ing reflected by a set or inteinal 
state variables (lSV‘s) wmen enter tne proolem 
oe scr i pi ion via the constitutive equations. 
Thus, tne necessity to moae'. each c-'ack -'to ter: 
{or nundreosj of finite elements is onviated with 
little Jos-- of accuroC) an; considerable 

computational Sowings. : u : m v* er.j, tne approach 
is sir;: \-r se If -cons i stent schemes to! o r -j 
CiObd’. i.-vw- luctr.OJS -t’. i'*C = C in Ctr.e: 

aOp • - u — . ■ uj < 

Tne necessary p s rt: c* tne rooel £’*e as 
follows: a Kinematic ascription c: tne aamage 

state; 2) c damage dependent set c: ply level 

stress-strain relations *.MCn account for matrix 


cracking; 3) 

a carnage oepencent lami nat' 1 on tnecry 

wn i c f ■ iiiuu- 

t s tne 

-erects of lhterply 

oe 1 ami na t i o; : 

; *) c se: 

■ carnage evolution lews 

fc r predict r 

•c me loe _ ’ 

»_to r y dependence c* tne 

carnage it ait 

a: eacr 

material point: 5} a 

Structure! j 

igoritnn fu? 

noaelling tne response 
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of components with spotlit 11.) variable and load 
history dependent stresses and damage; and 6) a 
failure function for predicting unstable crack 
growth leoding to loss of component structural 
integrity. Tnese parts are described briefly in 
the following subsections. 



c aooaoc 6D000C m —.": ;. 2 oooc^-pi 

LjkmJ CjtcIm <0 

Fig. f. Deiamination Area vs Load Cycles Grapr,: 
(0/9o ? ) c Specimen 





rig. 6. Normalized axial Modulus vs Del ami nation 
Area Graph: (u/ 9 G?) s Specimen 



seem tr.at tney should be modelled first. 
Typicaiij, these are o: a scale wnicn is very 
small compared to the geometry of a structural 
component. inus, it is hypothesized mat a local 
volume eienier.i may be seiectec which is small 
Compared tC the bGuriucfy vdue problem of 

interest an.:, at leas: for the case of matrix 
cracking, tr.~ damage can be assjmec t c be 
statist ice i i„. spatially homogeneous in mis 


element. Hoceiling the effects of matrix cracxs 
would then appear to be ideally suited to the 
approach taken in continuum damage mechanics. lr. 
this approach, first proposed by L.K. Kachanov in 
1958 |8). it is hypothesized that the effect of 
microcracks he locolly averaged on a scale 

whicn is saw 11 compared to the scale of the 
Structural component. Although tne procedure has 
been extensively uti slued in the literature, 
until I'ecci.t it r.a- i»wt beef: cbplitJ tC 

i oth ifjatrO u i '.hut TOP 1 v. i'>e J • - i - • * J 1 • 


over oC ! Hi t hr *■ I neHic t i ■. v ! ! t v t i C - ' c •- *• - '■ - 1 • 
me luCc i vi Jiic c 'cnic’:! ..OS ta*fc" Pj ,fl ’ u Ti" . 
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M 

where „ . ss me 1SV for matrix cracking in each 
ply, V L is on arbitrc-'i ly chosen local volume 
element or p ly thickness which is sufficiently 

large that a? . ooes not depend on the size of V. , 
c ^ 

u. are crac* opening displacements in V^ t 

u ^ are tne components of a unit normal to the 
crack faces, ana £ c is tne surface area of matrix 
cracks if 1 V , q ^ srtuwii it-. ' iy. a. 



Damage Dependent laming i ion Theory 

Unlike the ply level model for matrix 
cracking, statistical homogeneity cannot be 
assumed for delamination;,. Although statistical 
homogene i ty appears to hold in the plane cf the 
laminate, me same cannot be said tor tne 
tnrouen- in ickness direction. lnjs. da's. age is 
accounted for *ia area a vr»cging in tne laminate 
plant, cccuiiv-jnicd l .• »ineme:i: ti sumpticr- 

inrujcr, ir— li iCkness. i ng !j , tne iaminjtc 

eace: t > ■:t>nsi r -:tr. i 





rotations =. u a no „ ^ for the i t r. oe : ami naiec ply 


interface, a: snown ir. rig. 9. Kathemat 1 ca Uy , 
then ] 1 7 i 


fig. £. her ormet 1 on GeOfitctry for Region A^ 


P )s 




i ors 


f - o r 1 1 . « L.cse w ne r c cii matrix cracks arc 
normal ic t n- laminate pi-ne the local volume 
averaging process can oc shown to resuit in tne 
following p i ) level stress-strain relations 
{15,16!: 



M 

u 

yy 

r 


y.hc r c . o?'c component i or tne etestic 
1 J 

(undamaged) modulus tensor ir. ply coordinates, 
and tne subscripts L imp!., tiioy the components of 
me stress arc strain ter.scr are locally averaged 


y t 

• u ih>) - - it - 


1 - *<:-:*> 
(3) 

• V 

* > i*.y) - - 1* * 


* «(*-:,) S 
(<) 

ins 



• w"(x,y) < 

X 

(51 

wner * the supe^scn pt s 


imply undamagec 

micsurfcce quantities, 

end 

H ( 2 - : , ) is tne 


Heavy-sice step function. Also, a repeated index 
i in a proouct is intended to imply sunn.at ior., 
and the superscripts 0 imply displacement 
components across the dr iami nation. 



fig. 9. Schematic c-f ue laminated Region in a 
Composite Lay up. 
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The g lacement equations are averaged over 
the 1 oc a ) oreo in oroer to produce locally 
.averaged displacements tc be utilized in the 
I laminate formulation. Thus, 

vi L (».y.<> •• i ! Iu‘- Jit*’ H(I-I ,>(*“) ) * u'ldidy 

1 (6) 

,(«.*. i) * i I It'- »l»‘. «(*-*,>(*,) ) * «<J-J,» 

‘ * L (?) 

*n<3 

- ; J I*." - HI:-;,} 1 o *d> (S) 


b> averaging tne a ) sp i acenients , me ae lamina t ion 
jump discontinuities dre also averaged 
over A . 

Tne rc$u '• t ant laminate equations may be 
ootained d, suDstituting eolations ( 6 ) through 
( 6 ) into tne stra in-q ij>p iacement equations, end 
mis result into eq«a: i»-n: [:). T r. i s result i: 
tneh iniegiotcd tnrouu> me laminate thickness 
ano me divergence metres .s employed to oDtcin 
tne t o T sowing laminate editions 1171 : 


Ki - jjlOI, C k - ..f> • 1 j, l51 k * : i-l> u t> 


} 

1 »A 


. >‘ C V, -i / 


f o' 
i 0 

c*i _ j 0 

.1, f 2 i ’ 2 M* ^ 2 * i / D 


c 2i 
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4 i 
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“31 



L c J 
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5 \'*U (2 k * \-i> - \ 


(5) 
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l .» ! 

^°J 
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where |Hj and [H] are the resultant forces and 
moments per unit length, respectively, 

and io M l k and Uj) represent the damage due to 
matrix cracking and interply delamination, 

respectively. furthermore; , n is the number of 
plies, and d is the numner of delaminated ply 
interf aces. 


T he 


internal 
0 


state 


variable 


f or 


de lami nat ior, , u^ ; , is obiained Dj employing the 


divergence theorem or c 
me laminate. 

Tne rest It me proceo-’ r 


i o c a 1 vo 1 ume element e* 


J 2 i 


b 

" 2 -. 


I i 


v°r. d£ 
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01 
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, . n ob 
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(lie) 
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01 A. x 1 Z 
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where tne subscript i is associated with the ith 
delaminated ply interface. Furthermore, is 
equivalent to t^A^, where t^ is the thickness of 
the two plies above ano oelow the aelami nation 
1 17 | . By aef inition, the 2 component of the unit 
normal, n.,, is equivalent to unity. The matrices 
JQl with subscripts k are tne standard elastic 
property matrices for the uncamaced plies. The 
matrices !Q| with subscripts i apply tc the ith 
C6 i ami r.a tcO p • y inter* ace. snej represent 
average prourrties of tne plies above ana beiow 
the de lamina; ior. 117;. 

While tne oama ge -u cperiuc n i laminate ane'ysis 
moael may oe used to predict any or the effective 
engineering moduli for a laminate, experimental 
results are only availctie rcr tne axiai moduiLi 
and Poissor. ! s ratio. Therefore, the genera t 
utility of me model is demonstrated by c om o a *■ i nc 
model p^eoici ions to experimental results ror z x 
ano 


x ; 

Moo el 




been maoe 
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typical graph ite/epoxy system. The bar chart 
shown in Fig. 10 compares the model predictions 
to the experimental values for the engineering 
modulus, F , for combined matrix cracking and 
delamina! io:; the deldminrt*. ion interface 

location and percent of delamination area are 
listed in me figure underneath the laminate 
stacking sequence. As can be seen, the 
comparison between model results and the 

experimental results is quite good. Some limned 
results for ^0 is son's rj; i c d^e giver In ric. 1! 
using tnt Scmr 06 ** Chert fcrmat. hi;r, tr •_ 
exception or me 10/9G-, j,. laminate, tnese results 
are also quite good. Note that these results 

ha vc been obtained for c single set of input oata 
whicn do not depend on stacking sequence. 

Evaluation of Plv Stresses 


In crae** 

to 

evaluate tne stress 

state in 

each p i , it 

is 

firs: necessary to 

substitute 


OiTiaje Deoer>CBni Laminate E, 



bei. interlace 0/9C 0/90 0/9C 0/90 -*5/S0 *4S/*45 

expennerta! riooe: prediCUOft 


Fig. 10. Comparison or Experimental Results and 
Mcuel Predictions of the Laminate 
Engineering Modulus, F y , Degraded by 
bo: n Karri x Cracking ano Del ami nation 
i; ariiaq e . 


di sp laceweni equations (6) through (8) into the 
locally averaged strain-displacement equations. 
Utili 2 ing the divergence theorem on this result 
will then give the following equations for the 
strains in each ply 1 19 | : 


£ l * L i - 7 

A A 

U L ♦ ) t.^| - H(2-2 . ) 
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D*p#fxS*fi1 Lamnale v iy 

ix_ 



[0/901, |C. /9C ? |, »C 1Z 3 ) g tC^;, 

% D»(«mm*tior, 1€.6X 49. SV 24.2% 3S.3N 

Da:. tm«r1«ea 0/60 0/90 0/90 C/9C 


Fig. 11. Comparison of Experimental Results and 
Kooel Predictions of the Laminate 
Engineering Poisson's Patio, ■. J{v , 
Degraded b.» botr. Matrix cracking and 
Dei aminat icr. Damage. 
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The above equations »ay be utilized to 
obtain the ply strains, and these results may be 
substituted into equaiion (?) to obtain the 
stresses in each ply 1 19 j . Since the ply 
stresses determined by this procedure represent 
locally averaged values, they must be considered 
to be far-field stresses, so that stress 
intensity factors would ue needed in order to 
determine matrix crack-tip stresses. This point 
will oe discussed 'urine:* in the section or 
damage evolution ld*s. 


»ust be considered to be far-fleld stresses, so 
thata 1 SV evolution lows are generally of the 
form 1 19 ): 


Mj * °ij 


<l ki’ T ’ 
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°kk ' 
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'1 1 1 J 


(18) 


anu 


.0 
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i j Vl * 



(15) 


Z- COfl-hulc' COdr r.c: brr ^ COnSl ruC’.t." 

determine tne effect o: Jomage on the ' far- rieic* 
ply stresses in composite laminates !!?!. 

Predicted stresses are snown ir. Figs. 12 anc 
13 for a typical crossply laminate and a 
candidate guas i- isotropic layup with representa- 
tive damage states. *•$ evioenced from the 
results, me damage s igni f icent 1> affects the 
far-fieid ply stresses. Matrix cracks had a 
significant effect or. ply stresses in the S0 : 
plies in me cross-p';. laminate. The auasi- 
isotropic laminate exhibited a small stress 
reduction m the z4$' plies, but showed a 

substantial stress reouction (fifteen percent 
verses one percent) in ;ne 90“ plies. 

Dameoe Evolution Laws 

Since me ply stresses determined by this 
procedure reoresent local averaged values, they 


„ne r c * • , * , ano n i * ■ are t ne stress intensity 

fflClori, relate me ic r -rietc btr^sses to 

the crack t ip stresses ter a giver, c ■" a c ► 
geometry. however, it is assumed mat tne 
qeometry of both matrix cracking and 

oe lami nations is sufficiently independent ot 
stacking stouence tnat me stress intensity 
f actors may i.e treated as “'materiel properties” 
anc thus possess tne same stress intensity factor 
dependence for all steering sequences* inus, 
they are encompassd implicitly in tne material 
constants reouired to characterize damage 
evolution laws (IS) anc (19). 

One approach to me formulation of the 
internal state variable evolutionary 
relationships is through micromecnanical 

considerations. However, mis approacn is 
dependent on tne availability cf micromechanical 
solutions tnat can mooei tne essential physical 
character! st ics of the damage state. For the 
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problem of matrix cracks embedded in an 
orthotropic medium that is layered between two 
other orthotropic medic, the solutions that are 
currently available are applicable only to very 
specific loading conditions and damage 
geometries. Inerefore, the evolutionary equation 
proposed herein is phenomenological in nature. 
The form of inc damage r vw lul iunary relationship 
employed i r- this paper is based on tne 

Observat ior. made Dy h'or.g, cl cl. !?0| that fo r 

some iiidu-i.r. :s tne rate O' carnage surface 

c ^ 

evolu: 'c r i r ■ ! cac stes. — ~ , follows a ro«e*" 

ia* in me enr!'s_. release rate, C. Thus 1 2 : ! . 


amplitude of 26 ksi is sno*n' in Fig. K, 

Good agreement is found between the mocel 
predictions and the experimental results. 

To examine the amount of stress 
redi str ’but .or. that occurs during the damage 
accumulation, the model was used to determine the 
axial stress m the 90 plies c f tne ! 0^/90^ ! ^ 
laminate ioi*guc loaded at three cirrereni stress 
ampl iiuccs. * ig. It vruws tnat ter tne stress 
ampiUuue c: Jt ksi, the <1*10’. stress in tne 90' 
pi-cs a’ it*! Tor'*. j inuusd?. c.c*es j e$s mar. 

• ' . * ; n* or St m : r » c :■ r * r-v 


-• d -Vi 



K 

The term d^.^/dS reflects the changes tc 

tne interne' state -.a-*] able with respect to 

changes 1 ^ the damage surfaces. dci^/dS can be 

ootainec cnclyticall^ from relationships 

describing ;:. c kinematics of the crack surfaces 

for given damage states and loading conditions, 

snould Sod:, solutions exist. For transverse 

matrix creos in cress;:, laminates, tne average 

crack face displacement in me pure opening moae 

can be appro Aima ted Oj c solution obtained by 

Lee, et a.. \22\ for c medium containing an 

infinite number of alternating 0' and 90” 

plies. Thus, da.j/db can be determined for 

crossply laminates subjected to uniaxial loading 

conditions. It nas been found mat for typical 

continuous fiber reinforced graph ite/epoxy 
K 

systems 0 ^. ./dS can oe assumed to be constant 
for a given applied load until the damage state 
has reached an advanced stage of development. 
This assumption has facilitated the determination 
of tne material parameters • ^ ano n. Tne 

material properties *o r 4^4/3501-6 graphite/epoxy 
are useo in tne calculations to enable tne 
comparison of model prediction tc experimental 
measurements made by Cnou, et el. J24). The 
material parameters fc ! * this polymeric composite 
system nave been found to be 



Fig. 14. Katrix c rack damage in the 90" plies 
of a [0 ? /9G-J, AS 4/3501 - 6 laminate 
loaoed at a stress amplitude of 26 
s i . 


•i.hU c.n.; 


r = 6.35 (21) 



The damage history for a typical crossply layup 
has been predicted using 




prediction re* the oam=ge state ir. tne JC->/5ui; s 
laminate fatigue ioaoec at a maximum stress 


ucmage induced stress redi stribut ion 
*. r. t n p Of. pile 

1 cri nates subjected ' 
cmolitudc feticue loadinc. 


r, C * 3 " * C ' ^ C. C r 
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undamaged laminate. These results demonstrate 
that the stress redi stri but ion characteristics 
among the plies in the laminate are dependent on 
the loading conditions. These redistribution 
characteristics will affect the manner in which 
damage develops in the surrounding plies as well 
as eventual failure of me laminate. 

life P red ict ion 


redistribution as damage develops for a given 
fatigue luoJ history. 1ms predicteec stress 
redistribution in turn affects tne evolution of 
damage, thus producing a life prediction model 
which cen pc used for any stacking sequence 
regardless of the load history applied to the 
component. However, while ini t id moo el 
comparisons to experinicf.i are favorable, further 
research u suggested before the mood is 
ul • 1 5 .'rj i r. *. be s '< gn b c l . ■ r -s • 


Usually, ultimate ‘ci lure c f laminate: 
composites ;s cause: r., large scale fracture 
which is i nausea i-j fiber fracture at 
ae lamination sites. Therefore, tnere is not only 
a synergistic effect between matrix cracking and 
de lami nat i on , out also between delamination and 
me ultimate failure event. Since fiDer fracture 
occurs very near tne er.j of tne component life, 
ratner mar moael it with an additional internal 
variable, is is preferrobie to simply treat it as 


tne u 1 t i 

in»atc fail 

ure e^ent a no moo e 1 i 

L witn a 

fa i lure 

f u nc t i On 

. Typically, one 

Cf uc 

approach 

es C 0 u 1 ci 

be taxer. A pnenome 

no logical 

D i C c l. n 

-ou Id entail tne rood if icatic 

in of ar. 

existi nc 

f a i 1 ure 

function suen as the 

Tsai-Wu 

f a i 1 u r e 

f u nc L i on 

to account for the 

existing 

damage 

state. 

Alternatively, a 

f racture 


criterion could also br fioaified to account for 
tne carnage induced stress red istribut ion. The 
authors are pursuing mis subject further at this 
time. 

Cone 1 l- s i on 

The authors nave sno*n met Dy construct inc 
local overages of the r. inematic effects cf 
m l cr oc r a c x i r< g it is possible to construct 
continuous ' internal variables which appear 
explicitly in a modified laminations theory fer 
layered composites witn demage. Comparisons of 
predicted stiffness less as c function of damage 
state to experimental results lend creoence to 
tne moael. 


secduse me lamination tneory is damage 
aeoenoent . it produces predicted stress 
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APPENDIX C 




LIFE PREDICTION IN LAMINATED COMPOSITES 
USING CONTINUUM DAMAGE MECHANICS 


by 


David H. Allen 

Aerospace Engineering Department 
Texas A&M University 
College Station, TX. 77843-3141 


INTRODUCTION 


Laminated continuous fiber composite structural components are observed 
to develop substantial load induced noncatastrophic microstructural damage. 
Since this damage quite often leads to very small changes in linear elastic 
material properties, analysts are sometimes lead to the fallacious assumption 
that linear elastic analysis need only be supplemented by a knock-down safety 
factor in order to design away from failure. This supposition is erroneous. 
Although substantial globally stable microcracking can be accepted without 
component failure, the progression of this damage ultimately determines the 
life of the structure. Furthermore, the damage is spatially variable, tending 
to congregate in regions of high stress, and since stress concentrations are 
unavoidable in heterogeneous media, the development of damage is 
inescapable. In fact, this process can be considered to be in some sense 
beneficial because the integration of two separate brittle materials results 
in this damage, which is a ductile-like phenomenon, thus causing load 
transferral not unlike that which occurs in metals. 

A complicating feature of life prediction is the fact that the mode of 
failure in a component of given structural shape is observed to be dependent 
on the load history involved. In the same way that monotonic tension and 
compression in composites is observed to result in differing failure 
characteristics, varying the fatigue spectrum applied to a particular 
component car, result in gross changes in the observed failure characteri st ics. 
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Most current methodologies for life prediction in composites seem to be 
extremely limited in nature in that they are largely ad hoc in nature; that 
is, the nit-thods of analysis are in most cases dependent on component geometry, 
load history, and even stacking sequence and environmental cor. ' i t ions. 
Probably the most direct procedures utilize phenomenological failure laws not 
unlike Miner s rule | 1 , L' ! . As a result, they apply only to predetermined 
geometry and stacking sequences. 

Perhaps the most ambit ^us approach attempts to bring all of the details 
of the mechanics into rocus by using linear elastic finite elements and 
treating eacn crack as an internal boundary subject to growth determined by 
fracture mechanics criteria |3,4|. While this is desirable and probably the 
most accurate approach, it would require supercomputing capability for a 
typically complicated damage state such as that found in the circular cutout 
shown in Fig . 1 . 

Ine current chapter of this text reviews an approach to this problem 
which, although not complete at the time of this writing, has the potential to 
model the life of a laminated structural component, no matter what the 
geometric shape, and for any loading history, given only a set of input aata 
which do not oepend on the stacking sequence. Furthermore, the structure is 
modelled as a simply connected domain, with the effects of microcracking 
reflected d> a set of internal state variables ( I S V ‘ s ) which enter the problem 
description via the constitutive equations. Thus, the necessity to model each 
crack with tens (or hundreds) of finite elements is obviated with little loss 
of accuracy a no cons iderau !e computational savings. To this end, the approach 
is similar to self-consistent schemes [5| and global-local methods |6| 
utilized ir, other applications. 

Although the author is not so bold as to suggest that the methodology 
proposed herein represents a panacea, it is hoped that this approach is at 
least moving towards a generalized analysis tool applicable to laminated 
composite structural component life prediction. In the same way that 
classical plasticity theory developed in the early part of this century for 
isotropic and homogeneous media, the model possesses some features which are 
already grounded in fundamental mechanics, and some parts which are 
phenomenological at this time and will thus require future refinement. 
Therefore, the model should be treated as in a transitory state, and each user 




will want to supplement model predictions with a liberal experimental 
verification program. The usefulness of the model will thus vary from one 
application to another, am will be measured by the number of experiments 
which can be avoided with model usage. 

1 he primary objective of the damage model is to develop the capability to 
predict the fatigue life of laminated structural components. Accordingly, it 
is essential tnat the constitutive inputs to the model be independent of me 
geometry ano stacking sequence for tne particular application at hand. The 
approach taker, to this objective is to construct a nonlinear damage dependent 
lamination theory which can be implemented to any computational structural 
algorithm such as a finite element code. This then allows for modification of 
linear elastic codes via a time stepping scheme to account for load history 
dependent damage and an iterative algorithm to account for nonlinearity on 
each Lime step. The necessary parts of the model are as follows: 1) a 
kinematic description of the damage state; 2) a damage dependent set of ply- 
level stress- sc ra in relations which account for matrix cracking; 3) a damage 
dependent lamination theory which models the effects of interply 

de laminations; 4) a set of damage evolution laws for predicting the load 
history dependence of the damage state at each material point; 5) a structural 
algorithm for modelling the response of components with spatially variable and 
load history dependent stresses and damage; and 6) a failure function for 
predicting unstable crcCx growth leading to loss of component structural 
integrity. The chapter i-.tii review recent developments in each of the above 
areas, and support for tue model will be demonstrated via experimental 
results. 

MODEL DEVELOPMENT 

As shown in Fig. 1, a typical laminated continuous fiber composite with 
brittle matrix usually undergoes three significant and physically different 
modes of damage: matrix cracking; interply delamination; and fiber breakage. 
Normally, a substantia! distribution of matrix cracks occurs before 
delamination initiation, and the culmination of this phase of damage is 
accordingly termed the characteristic damage state 1 7 ) because it appears to 
reach this configuration independently of the load history. Del aminat ions 
then begin to appear as a result of stress concentrations at the i nter sect ions 
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of matrix cracks in adjacent plies, as shown in Fig. 2. Finally, the 
development of large scale del ami nations causes stresses sufficient to induce 
fiber breakage, and this usually leads to component failure via unstable 
fracture. A significant literature exists detailing the experimental 

observation .f these events | 8 - 1 2 ! . 

Sinct. iii.o. i- cracks appear to be the first mode of damage to occur in all 
cases, u aou.-.i seem that they should t»e modelled first. Typically, these are 
of a scale, .i ich is vei _> small compared to the geometry of a structural 
component. .• example, typical Saturation crack spacing in graph i tc/epo 

is about : irt, cracks per inch (for a ply thickness of 0.0055 in). Thus, it 

is hypetnesiicd that a local volume element may be selected which is small 

compared to : tie boundary value problem of interest and, at least for the case 
of matriA cracking, the damage can be assumed to be statistically spatially 
homogeneous in this element, as shown in Fig. 3. Modelling the effects of 
matrix crack v ,\Ould then appear to be ideally suited to the approach taken in 
continuum damage mechanics. In this approach, first proposed by L.M. Kachanov 
in 1958 | 13). it is hypotnesized that the effect of microcracks may be locally 
averaged on scale which is small compared to the scale of the structural 

component. - I though the procedure has been extensively utilized in the 

literature, until recently it has not been applied to laminated orthotropic 
media | Id- 18 , . 

Kinematic Description of the Damage State 

A st raignt forward and direct approach to averaging the kinematic effects 
of cracks within the local volume element was taken by Vakulenko and M.L. 
Kachanov in lb/1 1 191. This average is given by the following second order 
tensor: 


M 


u 


i j 



c c . 
Vj dS 


( 1 ) 


where J.' . is the 1SV ur matrix cracking in each ply, V^ is an arbitrarily 

chosen local volume element of ply thickness which is sufficiently large 

M c 

that u-j does not depend on the size of V^, u. are crack opening 
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displacement*, in V^, n^ are the components of a unit normal to the crack 
faces, aru! i is tne surface area of matrix cracks in V^, as shown in Fig. 

A. It t t > De shown in the next section that definition (1) is 
thcrmodyr.ami. .1 i ly consistent with the constitutive equations developed herein. 


Ply Lev el tit re as - Strain relati ons 


(Jons '.dr.-, '.ne local .< iuiiu: element shown in Fig. 5 with traction boundary 
conditions the external surface Sp In addition, the interior is 
assumed to Pu composed entirely of linear elastic material and li.'.w ich 
may include Lhin surface layers of damage). Integrating the pointwise 
Helmholtz free energy over the local volume will result in 1 20 1 


hr, - A, -t B . . .<.. • . .. . . c i , ,+D, aT. + E. . .c. . .aT. -+ i 2 F aT (2) 

EL L l t j Lij iijkl LlJ Lkl L L LijLij L L L v ' 


where A^ , 6,^, C L ^j k1 , u L , E j , and F L are locally averaged undamaged 

material constants. Furthermore, the locally averaged stress is defined by 


'Lij 


V 



and the loca averaged strain is defined by 


(3) 


«- Lj j i i ■ A a(u -n . t u j n i )dS (4) 

' S 1 

where n^ are components of the unit outer normal vector to the surface Sj. It 
follows that the governing field equations may be integrated over the local 
volume to op: a ■ n 


1) conservation of linear momentum: 


o, . . 

L J i , J 


( 5 ) 
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2) conservation of angular momentum: 


°Uj * Lji 


( 6 ) 


3) conservation of enery, 


u t - ' i ; j t- i j + q Lj, 


(7) 


where lj ‘ is as defined q^j is the locally averaged heat fij>, vecto* 


is the locally averaged heat supply; and 

4) the entrop> production inequality: 

r 


>L “ T 


* 0 


( 3 ) 


where s^ is trie locally averaged entropy, and is the locally averaged 

temperature. The effective locally averaged internal energy, u^, is given by 


U L U r.! 1 U [ 


(9) 


where u^ represents ths 
uncracked body, given by 


.ally averaged internal energy of the equivalent 


• c 

U EL = 

\. 

udV - 


i- 

V 

L 

where it 

a f r 

Called 

uncracked 

pi . 

acting 

mechanical 


r output 


v l s ; 


r i u i dS 


( 10 ) 







(ii) 


where 7 L arc fictitious tractions applied to the crack faces whicn represent 

1 £ 
the difference between tn-r actual crack face tractions and . Equations (5) 


ORIGINAL PAGE IS 
OF POOR QUALITY 




through (8) arc identical in form to pointwise field equations in a continuum 

I20|. 

On the oasis of this similarity the locally averaged Helmholtz free 
energy is now defined to ne |2G|: 


where it can he seen from equation (2) that h^ is the locally averaged 
undamageo elastic tic ..tw.. itz free energy. 

1 he similarity between the pointwise and local field equations leads to 
the conclusion that |21| 


,.n 

1 j a 

Using equations 


j 

( 1 ) 


oh 


El 


b i 


Li J 



(13) 


and (11), along with Cauchy's boundary conditions gives 


• c 
U L 



(14) 


where ,;f.. 3 :'e local Ij averaged equivalent crack face stresses [ 20 j . Thus, 
the thermodynamic significance of definition (1) is established. 

Now consider equation (14) in further detail. The kinetic quantities 
of., may he interpreted as generalized stresses which are energy conjugates 
to the kinematic strain-like internal state variables . , • It can be 
inferred from inis that tnere exists a constitutive relation between these 
variables or" tne form 


c 

Li j 


kl 



(15) 


which depends on loading ni story via the explicit inclusion of the internal 
state variauie-. . 

Therefore, substituting (15) into (14) will give 




■ " ] • c 

J uf (t )d i 


u l. (l Lkl^l 


\(ti). 


I' 1 Lkl 


(tj)) 


( 16 ) 
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Thus, expanding equation (io) in a Taylor series, substituting (2) and (16) 
into equation (13) and neglecting higher order terms yields |20|: 


1 M 


l ■ j 


"[ i jk l l Lk 1 i jk l u kl 


(17) 


Restrict iny nr. uamage to small quantities constitutes a sufficient but not a 
necessary coni it 'ion for dropping tha higher order terms. Equations (17) injj 
be written in t.ne following incremental form for isothermal conditions 


do... - i,' ., ,di, , 

L l j t.tjkl kl 


(18) 


where 


'Li jk 1 


_LlJ. 
' " Lk 1 


(19) 


Tnus, it can be seen that the constitutive equations can be represented 
equivalently by both equations (17) and (18). In the first case, the internal 


state is represented by a strain-like variable, a.j, whereas in the second 
case the effect of damage is incorporated within the stiffness tensor, 
C i j k 1 ’ I s implicitly damage dependent. A careful inspection of 

equations (17) and (18) will show that 


C!,, can be determined as a function 
i Jk 1 

by differentiating (17) with respect to 


of the damage internal v,-.r iobles a. 

I J 

the stress tensor. :rms, it is found that C 1 j ^ ^ is a function of 

which is a fourth order tensor. Obviously, either 

could be termed the damage variable. In this work, tne 

author prefers to call { ■ the damage internal variable, and tne 

measure of the damage state. This is due to the fact that in brittle systems, 
it can be seen from equations (1) that when the damage surface area is 

3a ij /jC Lkl is 


3a . ./3 l , ■ , 

TJ LlJ 

a . - Or 3;. . ; 
1J 1 J 


, changes linearly with strain, so that 


constant , 

constant. Thus, the latter is a better measure of the damage state. 

Equation (17) and tne associated internal energy have previously been 
shown to be exact for both a linear elastic isotropic cracked body |22| and a 
linear elastic orthotropic cracked body 1 23 1 when the cracks are evenly 
spaced. For the case where all matrix cracks are normal to the laminate plane 
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coordinates. When the matrix cracks are not normal to the midplane, a more 
complex procedure is required [25] . Fortunately, this circumstance appears to 
arise only wnen several adjacent plies are stacked at the same orientation 
1251. 


Damage Dependent L aminat kt, 1 heory 


Unlike : ply 1 e * r i model for matrix cracking, statistical homogeneity 

cannot be assumed for de laminations. Although the delaminations appear to be 

evenly distributed in the plane of the laminate , the same cannot be said for 

the through-: nickness direction. The damage is therefore accounted for via 

area averaging in the laminate plane, accompanied by a kinematic assumption 

through the thickness. Accordingly, the laminate equations are constructed by- 

assuming tnat the Ki rchnof f-love hypothesis may be modified to include the 

effects of jump displacements u?, v*?, and wP, as well as jump 

L) D iii 

rotations a- and u,. for the ith delaminated ply interface, as shown in Fig. 

6. Mathemat ical ly , then 1261 


u ( < , p , - i u (x,y) 
v(x,y,z) - '/(x,y) 


+ H(z-z.) 

D I 
D i 1 

+ H (z-z i ) u? 

(21) 

+ H(z-z i ) 

,D, 
c, i 1 

+ H ( 2 - 2 i ) v? 

(22) 


and 
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( 23 ) 


o D 

w(x,y,z) - W (x,y) + H(z-2.) w. 

where the superscr ipts "o' imply undamaged midsurface quantities, and H(z-Zj) 
is the hca^iide step function. Also, a repeated index i in a product is 

intended to imply summation, and the superscripts D imply displacement 
components .1 :.ross the de lamination. 

T iv <■ 1 lac.ement equations are averaged over the local area, A , 
shown in 1 i.,. , in Greet to produce locally averaged displacements to be 

utilized 1 n : lie laminate \ ormu lation. Thus, 


u. (x.y.z) - x J |u 
L 


l P 


H ( 2 - 2 i ) (6^) ) + H(z-z i ) u*?|dxdy (24) 


V L (•’,.> 


1 v 


1 A, 


H(z-Z . ) (*•*?) ) + H(z-Z i ) v^jdxdy (25) 


and 


w L ( x . 


i . 


L A 


(w t H(z-2|) w- Idxdy 


(26) 


By averaging the displacements, the delamination jump discontinuities are also 
averaged over A . 

The laminate strains are given by 
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The i diii 1 1 ia ie constitution is obtained by integrating the stress in each 
lamina over r lit- laminate thickness. The resultant midplane forces and moments 
per unii i-.i.iii. u» region , in the laminate are thus given b. 



(33) 


and 


M 

x 

M 

y 

M 

xy 



IVJ 


(34) 


where t is t he laminate thickness. 

The resuitant laminate equations may be obtained by substituting 
equations (2-1) through (26) into equations (27) through (32), and this result 
into equations (20). This result is then substituted into equations (33) and 
(34), and the divergence theorem is employed to obtain the following laminate 
equations |26j : 
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where [Nj ana !M} are ir.e resultant forces and moments per unit length, 

M D 

respectively, and in j. mid iu.) represent the damage due to matrix cracking 
and interply ..ie lamination, respectively. Furthermore, n is the number of 

plies, and ,j i the number of delaminated ply interfaces, as shown in Fig. b. 

The internal state variable for delamination, lu^i, is obtained by 
employing tno divergence tneorem on a local volume element of the laminate. 
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The resulting procedure gives |26] 


D 2 - D 

“ii 1 r: “w dS 

! ’ 


l 2i ’ v' 


, I s 2i 


• - D n - 

v . M e. j 


(37a) 


( 3 ■ 


U 

a 3i 


i u .n dl> 
d 2 i 


(3\) 


D 

u 4i 


]_ 

A 

l 


D HC 

!}■ . n dS 

t 2 


2i 


D 

a 5i 


1 

A. 


D , c 
- n d S 

t 2 


r i 


( 3 7 d ) 


(37e) 


where the subscript i is associated with the ith delaminated ply interface. 
Furthermore, v ; j is equivalent to t^A^, where t^ is the thickness of the two 
plies above .nu below tiic del ami nation, as shown in Fig. S. By definition, 
the z component of the unit normal, n 7 , is equivalent to unity. The matrices 
|Q) with subscripts k are the standard elastic property matrices for the 
undamaged p.ies. The matrices [Q] with subscripts i apply to the itn 
delaminate.) r - f, interface. They represent average properties of the plies 
above and below the deiami nation. For example. 
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(40) 
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(41) 


where the Supt-i scripts A and B designate the properties of the ply immediately 
above and brio* the de lamination, respectively. These are described in 
further oeta • i in reference |26|. 


Determine! ion of Damage bependent Stiffness 


A measure of the accuracy of laminate equations (35) and (36) can oe 
obtained o, comparing predictions of damage dependent stiffness to 
experimental results. In order to do this, it is necessary to construct tne 
(stacking sequence independent) material parameters developed in the previous 
section. Ttu- loading direction engineering modulus, E x , and Poisson's ratio, 


xy 


of tne .aminate are defined by |27] 
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(42) 


(43) 


where t is the laminate thickness. 

For tne purpose of comparing the model predictions to experimental 
results, this development is confined to the case of a symmetric, balanced 
laminate with del ami nation sites symmetrically located with respect to the 
laminate midplane. For inis special case, the following expressions for t x 


ORIGINAL PAGE IS 
OF POOR QUALITY 






Implement at ion of equations (44) and (45) to predict the damage degraded 
laminate moduli requires the specification of the partial derivatives of the 
internal state variables with respect to strain for a given damage state. In 
the absence :i I SV evolution laws, the damage state must be determined 
experimentally. Expressions for the internal state variables have been 
previously developed |28j by employing energy principles. 

For example, in the case of matrix cracking in cross-ply laminates where 
only the opening mode or fracture is involved the following expression has 
been developed 



(45) 
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where m is tne number of consecutive 90" plies, p is the number of 0" plies, q 

is the number ..f 90“ plies, E is the initial undamaged modulus, and E 

x 0 1 
is the aamdwt- degraded modulus corresponding to matrix crack damage state 

. 1 he t «r rm in the parentheses has been determined experimentally from 

tesls or a graph i u-/epc»y laminate and is given by \27\ 

r 

-p 1 ! - U. 59569 - u. 061607 S t 0.04623 S 2 (60) 

x . i j 
u 

Other components are as determined in reference 127). 

While the damage-dependent laminate analysis model may oe used to predict 
any of the effective engineering moduli for a laminate, experimental results 

are only available for the axial modulus and Poisson's ratio. Therefore, the 

general util it, of the model is demonstrated by comparing model predictions to 
experimental results for and v for the fully developed damage states 
illustrated in Figs. 10-14. The delamination interface location has been 
determined experimentally and the delamination area has been estimated from 
the x-ray radiographs using an optical planimeter procedure. In both tne 

model analysis and data reduction, it is assumed that the delamination sites 
are symmetrically located about the laminate midplane and contain the same 

delamination surface area. 

Model predictions ha\e been made for a typical graphite/epoxy system witn 
properties shown in Tabie 1. The bar chart shown in Fig. 15 compares tne 

model predictions to the experimental values for the engineering modulus, E , 
for combined matrix cracking and delamination. The delamination interface 
location and percent of delamination area are listed in the figure underneath 
the laminate stacking sequence. As can be seen, the comparison between model 
results and the experimental results is quite good. Some limited results for 
Poisson's ratio are given in Fig. 16 using the same bar chart format. With 
the exception of the | G/ 9 O 2 J s laminate, these results are also quite good. 
Note that these results nave been obtained for a single set of input data 
...nich do not depend on stacking sequence. On the basis of these results it is 
concluded that the theory proposed in equations (35) and (36) is a candidate 
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I able 1. AS4/3502 Graphite/Epoxy Lamina Properties 

E n 21.0xl0 6 psi (144.8 GPa) 

Eo'> 1.39x10^ pi i (S.fS GPa) 

L. L 

V-..-) ( 1.312 

1 u 

G j 2 u.694x!0 6 psl (4.79 GPa) 
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for damaged laminates. 


Eval uat i on_ ot P ly Stresses 

In oide. 10 evaluate the stress state in each ply, it is first necessar., 
to substitute displacement equations (24) through (26) into the locally 
averaged strain definitions (27) through (32). Utilizing the divergence 


theorem on 

Hi i 

s resu 1 1 » I 

i 1 i then gi ve 

the following equations 

for the strains 

in each ply 

- 





l L “ 

X 

0 

1 L . 

- Z Ik - 

X 

“<‘- z 1 > “Si 1 

+ H(z-z.) u^. 

(51) 

1 L ( 

y 

t. , 

i 

i 

- z U L ’ 
y 

ili-’-.-j) .“,1 

+ H(z-z • ) 

(52) 

1 L “ 
2 

0 

'■ l . 

i- 

’ H(z-z i ) 

u 

“li 


(53) 

L yz 

u 
* L 

- • U L„, 

> “ J 

1 "( 2 - z 0 a 4 i 

' 

(54) 

£L xz = 

C 

‘ [_ 

, ' K.... 

,vZ xz 

, H(2-2.) 

1 

(55) 

c L = 


1 

1 — 
X 
<•< 



(56) 


The above equations may be utilized to obtain the ply strains, and these 
results may r.c substitutes into equation (20) to obtain the stresses in each 
ply 1 29 ! . Since the p;., stresses determined by this procedure represent 

locally averaged values, they must be considered to be far-field stresses, ->o 
that stress intensity factors would be needed in order to determine matrix 
crack-tip stresses. This point will be discussed further in the section on 
damage evolution laws. 

A computer code has been constructed to determine the effect of damage on 
the "far-fieiu" ply stresses in composite laminates 1 29 1 . Results presented 
are for a given laminate --train e = .01 (all other strains assumed to be 
zero). Damage variables v,ere calculated for matrix cracks in a saturated 
damage state using equations (38) and (39) assuming Ug = 0.0001 . The off- 
axis and 90 plies uu: the matrix crack damage terms of ug and u^, 

respective.'!/. Ko damage is assumed in the O' plies. Since the laminate is 
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subjected only to c xQ , is assumed to be the only delamination damage 
component. This term is calculated for an equivalent delamination area by the 
use of equation (40) witn = .00001". 

Hit I'irsu 1 is obtained from the model are shown in Table 2 and Figures 17- 
22 . As eviutM.ced by ti.v results, the damage significantly affects the far- 
-field ; st e-.ses. Dan.mjc variables were calculated b> using equations (38) 
through (40; to simulate damage existing in several previously tested 
laminates. Matrix cracks caused substantial stress reductions in ply stresses 
in the 50 piles in cross-ply laminates. For example, in the 50 plies of me 
I O 2 /90 2 i s laminate the matrix cracks resulted in a thirty-four percent far- 
field ply stress reduction. The two quasi-isotropic laminates developed 
different damage resulting in dissimilar far-field ply stresses. 1 ne 
[90/±45/01 s laminate exhibited little matrix cracking, thus producing only a 
small reduction in ply stress in both the 90" and ±45" plies. The |0/±45/90| s 
laminate exploited a similar stress reduction in the ±45" plies, but showed a 
substantial stress reduction (fifteen percent verses one percent) in the 90 
plies when compared to the 1 90/ ±4 5/0 ] s laminate. This alteration in ply- 
stresses snouid significantly affect the growth of new damage in tne 
composite. 


Damage Evolution Laws 


It is n /pc'tiies ized mat the evolution of damage in each ply is driven by 
the current stress state at tne crack tips within that ply [20). Tnis stress 
state depends upon the strains, temperature, and internal variables witnin 
that ply. In addition, ^ i nee these are locally averaged values, they must be 
regarded r.s far-field quantities, and the stress intensity factors should be 
included ir, order to account for the crack tip stresses. Thus, the internal 
state variable evolution laws should more properly be written in the following 
form: 


and 


.M .M , 

Vi ' °i j c Lk i 


kt’ u ksP I 


K t . K 


IP 


K m> 


.D .D , . M 0 „ Y „ x 

J i j (E LkP ' ’ ‘kt ’ u kt ’ K P K IP K IIP 


(57) 


(58) 
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where Kj, Kjj, and Kjjj are the stress intensity factors, which relate the 
far-field stresses to the crack tip stresses for a given crack geometry. 
However, it is assumes ihat the geometry of both matrix cracking and 
delaminations is sufficiently independent of stacking sequence that the stress 
intensity factors may be treated as "material properties" end thus possess l ne 
same stress intensity factor dependence for all stacking sequences. Inus. 
they are encompassd implicitly in the material constants required to 
character iiir damage evolution laws (5?) and (58). 

One approach to the formulation of the internal stele . aricbie 
evolutionary relationships is through micromechanical considerations. 
However, this approach dependent on the availability of micromechanical 
solutions that can model tne essential physical characteristics of the damage 
state. For the problem of matrix cracks embedded in an orthotropic medium 
that is layered between two other orthotropic media, the solutions that are 
currently available are applicable only to very specific loading conditions 
and damage geometries. iherefore, the evolutionary equation proposed herein 
is phenomenal og ical in nature. The form of the damage evolutionary 
relationship employed in in is paper is based on the observation made by Wang, 
et al. I 30 1 tnat for some materials the rate of damage surface evolution per 

• j ^ 

load step, follows a power law as shown below, in which the strain 

energy re'ieait- rate, G, ana a material parameter, n, serve as the basis ana 
exponent, respec t i ve ly . Thus, 


aS Dr n 
dN - PG 


(59) 


where P is a material constant. To develop internal state variable 
evolutionary equations in tne form of equation (59), . must be related to 


M 


i J 


q.., as defined by equation (1), 
1 J 


the surface area of damage. Since 
represents tne kinematics of the crack faces, the damage surface area alone 
will not be sufficient to describe the crack face displacements. Assuming that 
each crack in the material volume shares a common geometry and orientation, 
then the specification of the far field strains will complete this 
description. Tne rate of change of the internal state variable can therefore 
be expressed as follows |31|: 

M 

. dS 

( 60 ) 


d ‘1j 


dh 


du i j ^ 

dS dh 
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where the far field strains <are reflected in the term du^/dS which relates 
the changes .r, the internal state variable during the development of damage 
surfaces. Inns, using equations (59) and (60), the stable evolution of the 
internal stale variable is given by 


_ii 

dS 


k , b' ON 


( 61 '' 


T ne 1 1 


.lu M ./dS 
> J 


inflects the changes in the internal state variable 

M 

with respect to changes to the damage surfaces. dj.^/dS can be obtained 

analytically rrom relationships describing the kinematics of the crock 

surfaces for given damage states and loading conditions, should such solutions 

exist. For transverse matrix cracks in crossply laminates, the average crack 

face displacement in the pure opening mode can be approximated by a solution 

obtained oy lee, et al. 1 2 3 1 for a medium containing an infinite number of 

alternating (r and 90" plies. Thus, da.j/dS can be determined tor crossply 

laminates subjected to uniaxial loading conditions. It has been found that 

M 

for typical continuous fiber reinforced graphite/epoxy systems du.j/dS can 
be assumed to be constant for a given applied load until the damage state has 
reached an au.anced stage of development. This assumption has facilitated the 
determinat ion of the material parameters kj and n. The damage state at any 
point in the loading history can now be determined by the integration of 
equation (61) using the laminate averaged ply responses. This integration is 
performed numerically because of the nonlinearity of the damage evolutionary 
equation. The fourth order Runge-Kutta method has been found to be suitable 
for this app I ication. 


The development of t ne matrix crack damage state in crossply laminates 

subjected to uniaxial fatigue loading is examined using the proposed damage 

evolution equation. To maintain the thermodynamic admissibility of the 

fatigue loading process, it is assumed that the values of the internal state 

variables remain constant during the unloading portion of the load cycle. It 

is further assumed that trie required thermodynamic force is very small 

M 

compared to trie available thermodynamic force, thus a-- will change at the 
onset of load application. The material properties for AS 4/3501 -6 
graphi te/epex } are used in the calculations to enable the comparison of model 
prediction to experimental measurements made by Chou, et al. |32|. Tne 
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material parameters for this polymeric composite system have been found to be 
kj - 4.42 anc n = 6. 39 (62) 

Tne riamaqe histories for two crossply layups have been predicted using 

the mooe i . I tie model predict ions for the damage state in the I 

laminates fatigue loaded at maximum stress amplitudes of 38 ks- and 48 ksi are 

shown in \ igs. 23 and 24, respectively. 1 he lower stress amplitude is 

equivalent to eighty percent of the monotonic crack initiation strt-.s, while 

the higher stress amplitude is equal to ninety perent of the initiation 

stress. 1 he experimentally measured damage states were originally measured in 

M 

terms of the crack density. However, the corresponding for sach carnage 

level can be approximated bj the relationship proposed by Lee, et al. 123]. 

Good agreement is found between the model predictions and tne 

experimental results. The damage evolution for the thicker | O^/SO t ] s laminate 
is shown in fig. 25. This laminate was loaded at a maximum stress amplitude 

of 26 ksi. Ibis amplitude corresponds to eighty percent of the quasi-static 

matrix crack initiation stress. The results for this load case indicated good 
agreement with the experimental data. The effect of the load redistribution 
on the damage evolution is apparent in this load case. A markea decrease in 
the rate of damage evolution after fifty thousand load cycles was indicated by 
the model. On the other hand, the experimental data showed this decrease to 
occur after only ten thousand load cycles. Since the evolution of 
delaminations was not induced in the analysis, the predicted decrease in the 
rate of damage evolution is attributed only to the matrix crack induced 
transfer of load from the 90 plies to the adjacent 0'' plies and the resulting 
decrease in tne available thermodynamic force. The measured values of the 
damage state, however, may have been influenced by the formation of 
delaminations along the free edges and in the interior. Such occurrences can 
drastically affect the stress distribution among the plies and the available 
thermodynamic force for carnage evolution. 

To examine the amount of stress redistribution that occurs during the 
damage accumu iation, the model was used to determine the axial stress in the 
90" plies of tne 1 0^ /90-j i c laminate fatigue loaded at three different stress 
amplitudes, rig. 26 shows mat for the stress amplitude of 38 ksi, the axial 
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stress in the 90“ plies after forty thousand cycles was less than fifty 

percent of the original stress level in the undamaged laminate. Therefore, 

the rate of damage evolution is expected to be relatively low during the 

latter stages or the loading history. This is observed in Fig. 21 , which 

M 

shows the corresponding values of the internal state variable, u 7 p » for t ne 
three load cases. The 26 rsi stress amplitude load case, on the other hand, 
produced on!., gradual changes in the axial stress and damage state as o spared 
to the ctcu-r t stress amplitudes. The percentage decrease from the original 
undamage sire.,s level increased with the fatigue stress amplitude. These 
results demonstrate that tne stress redistribution characteristics among the 
plies in the laminate are dependent on the loading conditions. These 

redistribution character ist ics will affect the manner in which damage develops 
in the surrounding plies as well as eventual failure of the laminate. 


An Algorithm ; o r Structura : Analys i s 

In oroer to moael the response of a structural component with spatially 
variable stresses, it is necessary to incorporate the damage dependent 
lamination i. ivory into c structural analysis algorithm. This will be 
accomplished via the finite element method. To do this, first recall tne 
equilibrium equations for a plate [33]: 


jN 


3 a 7 y 


** - p 


(63) 


j N 3 N 

^ J 

a < hy 


a p 


(64) 


o M 


XV 


j x j y 


-0- = -P- 


(65) 


It is new advantageous to express equations (35) and (36) in the 
following niniru form: 
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( 69 ) 

(70) 

( 71 ) 

(72) 

(73) 
(79) 


Subst i lui mg equal mo (66) into equations (63) and (69) results in tne 
in-plane equilibrium equations [ 39 j . Similarly, substituting equation (67) 
into equation (65) results in the following governing differential equation 
for the out-of-plane deformations: 
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Integrating the governing differential equations, (75), (76), ana (77) 
against a variations in the displacement components and employing Green's 
Theorem twice results in the weak formulation for the laminated plate 
equilibrium equations. Inis form is described in detail in reference 34. 

The finite element formulation utilizes five degrees of freedom at each 
••nde, k. These consist of two in-plane displacements, u° and v°, one out-of 


plane displacement. 


and two rotational terms. 


x , y 
a. and of 


k’ 

The following 


k 5 — • J k k‘ 

displacement fields are assumed to represent the components of deformation 

within an element: 


u. 


e e 

i -U • 

= 1 J J 


J = i 


(76) 


0 
v , 
L. 


j = l 


e e 
v .v . 
J J 


(79) 


• e e 


(80) 


where i c , j - .w., 6 j , o : 
element, m is the number of nodes the element contains, and p is three times 
the number of nodes. 

Substitution of equations (78) through (80) into the weak formulation of 
the plate equilibrium equations results in the following system of algebraic 
equations for typical finite element: 


and represent the shape functions for the 
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where the sun ness and force matrices are as described in reference 34. 

The system of equations given in (81) can be used in conjunction with any 
type of element for which the shape functions are known. For the model 
developed in me present work, a three node triangular element with five 
degrees of freedom per node is selected [34). These in turn are assembled by 
the standard met hod into o global set of finite element equations. 

A fim it- element code has been developed based on the above 
formulation, ine implementation of tne damage evolutionary relationship shown 
in equation (61) requires that the solution algorithm be repeated for every 
load cycle. during each cycle, the ply stresses are calculated and used to 
determine i r.e increment in the matrix crack internal state variable for each 
ply. The updated damage state is then used in the calculation of the laminate 
and ply response at the next load cycle. 

The response of a laminated tapered beam subjected to uniaxial fatigue 
loading conditions has been examined. This beam has a 1 0/902 1 s stacking 
sequence ana possesses the material properties of continuously reinforced 
AS4/3502 graph i te/epoxy. Its length is 17.78 cm and the width is 4.50 cm at 
the clamped end and tapered to 2.87 cm at the end where the load is applied. 
The beam is loaded at a distributed load amplitude of 17.5x10^ N/m and R= 
0.1. Since the beam is symmetric about its length, it is sufficient to model 
half of the beam with a mesh containing 28 elements, 24 nodes, and 120 degrees 
of freedom. ine development of matrix cracks during the loading history in 
the 0“ and in.; plies, as represented by the internal state variables, is shown 
in Figs. 28 and 29, respect i vely . For clarity, the damage levels in each ply 
are normal izera by the largest value of the internal state variable within that 
ply at the end of 22495 load cycles. The predicted results indicate that the 
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matrix cracks first occur in the 90" plies at the narrow end and progress 
toward the wider end as loading continues. Axial splits in the 0° plies do 
not develop until after tne appearance of the transverse matrix cracks in the 
90' plies. iniiS, the transverse matrix cracks accumulate and loads carried by 

the 90 pi ie:. are transferred to the O' plies. Depending on the amplitude of 

the faiiyut: load, more man hair* of the load initially carried by l ne 

undamaged plies can m transferred to the 0" plies. The additional loads 
coupled mm me stress . n..entrat ions caused by the transverse matrix cracks 
create suiia.hle conditions for the growth of axial splits. The results su 
show that tne progression of transverse matrix cracks along the length of tne 
tapered beam decelerates near the midway point. The stresses in the region 
beyond the midway point are most likely to be insufficient in sustaining 

additional damage. Instead, the intensity of matrix damage increases at the 
narrow end during the letter portion of the loading history. This process 

will continue until the matrix cracks have either reached the saturation level 
or when the laminate fails. The model predictions are qualitatively supported 
by the experimental result shown in Fig. 30. 

Life Preoi ct ;ci > 

Usually, ultimate failure of laminated composites is caused by large 
scale fracture which is induced by fiber fracture at delamination sites. 
Therefore, tnere is not only a synergistic effect between matrix cracking and 
del ami riot ion, but also between delamination and the ultimate failure event. 
Since fiber fracture occurs very near the end of the component life, ratner 
than model it with an additional internal variable, it is preferrable to 

simply treat it as the ultimate failure event and model it with a failure 
function. Typically, one of two approaches could be taken. A 

phenomenological approacn would entail the use of an existing failure function 
such as the Isai-Wu (3b] and Whitney-Nuismer failure criteria 136] to account 
for the existing damage state. Since these failure criteria depend on 
stresses and si rains, respectively, in a given ply, they will automatically 
account for load history dependence, since the model predicts history 

dependence dur to damage evolution. Alternatively, a fracture criterion |37| 
could also be modified to account for the damage induced stress 

redistribution. In this case, the energy release rate would have to be 
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modified to account for me redistribution of stresses. The author is 
pursuing this subject further at this time. 


CONCLUSION 


I lie dull;..- nas shown mat by constructing local averages of the kinematic 
effects of mi e.rocrack ing it is possible to construct continuous internal 
variables wniui appear explicitly in a modified lamination theory for layered 
composites with damage. Comparisons of predicted stiffness loss as a function 
of damage statu to experimental results lend credence to the model. 

Because the lamination theory is damage dependent, it produces predicted 
stress reoistr ibut ion as damage develops for a given fatigue load history. 
This predicteeu stress red istribution in turn affects the evolution of damage, 
thus producing a life prediction model which can be used for any stacking 
sequence regardless of me load history applied to the component. However, 
while initial model comparisons to experiment are favorable, further researcn 
is suggested before the model is utilized in a design setting. 
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X-ray radiograph of a i 0 / 9 0 2 ]j AS4/3502 tapered beam specimen shov- 
ing the development of transverse matrix cracks. 


0 


ORIGINAL page IS 
OF POOR QUALITY 



APPENDIX D 




Damage Prediction in Laminated Composite? 
with Continuum Damage Mechanics 


David C. Lo* 
David H. Allen** 

and 

Kevin D. Buief 


Abstract 

A finite element model is developed for laminated composite plates ex- 
periencing fatigue load induced microstructural damage. Continuum damage 
mechanics based constitutive relationships, in which the formation of damage 
with a laminated composite plate is modeled with strain-like internal state 
variables, are utilized to obtain the weak form of the equilibrium equation. A 
damage evolutionary relationship for matrix cracking is incorporated into the 
model to provide the capability to examine the response history of laminated 
plate structures. Finally, this model is implemented into a finite element com- 
puter code and is used to predict the response of a simple laminated plate 
structure loaded in fatigue. 

Introduction 

The increasing use of laminated composite structures for engineering 
applications has lead to an increase in the number and types of methods in- 
troduced to analyze the mechanical response of this type of structure. Some 
theories incorporate the use of finite elements while others result in strictly 
closed form solutions. One of the earliest theories introduced for the analysis 
of laminated plates is known as the laminated plate theory. This theory is a 
modification of the classical plate theory that is based upon a set of assump- 
tions known as Kirchhoff's hypothesis. This theory is the simplest available; 
however, the out-of-plane components of strain in the plate are neglected. 
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This is due to the assumption that a line normal to the mid-plane prior to 
deformation will remain normal to the mid-plane after deformation. Mary 
higher-order laminated plate theories have been introduced in attempts to ac- 
count for the transverse components of strain. One of the simplest higher-order 
laminated plate theories is known as the shear deformable theory (Reissner 
1945 Mindlin 1951). Tne classical laminated plate theory is in fact a special 
case of the shear deformable theory. 

To this point, very little work has been performed on laminated plate 
theories that account for the damage formation within the laminate and the 
effect it has on the material response (Redd}’ 1987, 1990). Damage within a 
laminate first appears in the form of matrix cracks and for some cases can be 
the only type of damage to form. At the intersection of these matrix cracks 
with those in the adjacent plies, interply delaminations are initiated. Edge 
del animations can also be found along the free edges. These phenomena in 
turn contribute to the formation of fiber matrix debonds and fiber fractures 
which lead to catastrophic failure of the laminated composite. Experimental 
evidence indicates that the formation of matrix cracks and interply delami- 
nations have a significant effect on the material properties of graphite/ epoxv 
laminates (Norvell 1985, Georgiou 1986, Groves, et al. 1987). While other 
damage types affect tne ultimate life of continuous fiber laminated compos- 
ites, the formation of matrix cracks and interply delamination affect those 
properties that dictate the ultimate failure of the component. Thus, at the 
present time, it may be sufficient to consider only matrix cracking and inter- 
ply delarmnation in the formulation of a laminated composite plate damage 
model. 


New theories have been introduced in attempts to more accurately 
predict the response of laminated composites subjected to various loading 
conditions. One approach is based upon fracture mechanics; changes to the 
material behavior are characterized through alteration in local boundaries. 
Three dimensional finite elements are used to simulate observed matrix crack- 
ing patterns (Reddy, et al. 1987) and delamination formation (Murthy and 
Charms 1985, 1986). Since each crack is modeled separately, a large number 
of elements is necessary to construct an accurate finite element mesh. Mod- 
els of this magnitude will require a large amount of CPU time on a powerful 
computer, resulting in a very expensive solution to the problem. 

It is submitted that a more cost effective method for obtaining ac- 
curate predictions to the response of laminated composites is needed. One 
possibility is to use continuum damage mechanics. In this approach, the aver- 
aged effects of each type of damage in a small representative volume element 
are considered rather than individual flaws. Thus, the boundary value prob- 
lem is reduced to a simply connected domain with nonlinear material behavior 
and the number of finite elements is reduced to a manageable level. One of the 
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authors (Allen, et al. 1987a. b,c) has proposed one such continuum damage 
mechanics model and has produced accurate predictions to expei iment alh 
measured stiffness losses in graphite/epoxy laminated composites using the 
model. This model uses a set of second order tensor valued interna- state vari- 
ables to represent the matrix cracks and delaminations within the laminate. 
The incorporation of this theory into a finite element model will allow for the 
use of two dimensional elements in the modeling of a laminated composite 
plate. This will make it possible to model a laminated composite plate con- 
taining a large amount of damage with considerably fewer elements than those 
used by the previously discussed models. In this paper, the development of 
the finite dement model will be presented. The accumulation of matrix crack- 
ing in a laminated composite plate is then examined with this finite element 
model. 

Field Problem Formulation 

The formulation of the governing differential equation for a laminated 
composite plate with damage follows the same procedure as that used for an 
undamaged laminated composite plate (Agarwal and Broutman 1981). Dam- 
age in the composite plate is reflected through modified constitutive and dis- 
placement equations. The details of these equations as well as the develop- 
ment of the damage model can be found in several publications (Allen, et al. 
1987a, b,c, 1988). Therefore, only those portions of the damage model nec- 
essary for the problem formulation will be considered herein. The effects of 
the matrix cracks are introduced into the ply constitutive equations through 
a modification of the thermodynamic constraints placed on the field problem. 
Specifically, the Helmholtz free energy is modified to contain the mechanical 
effects of matrix crack formation. This form of the Helmholtz free energy 
is then used to obtain the constitutive relation for an orthotropic material. 
Symmetry constraints are then applied to obtain the following constitutive 
relation: 

K) = !<?]{-'! -an (1) 

where {& l } are the locally averaged components of stress, IQ] is the trans- 
formed anisotropic stiffness matrix for the ply, { sl } are the locally averaged 
components of strain, and {a™} are the components of the internal state vari- 
able for matrix cracking in terms of laminate coordinates. They are defined 
to be the matrix crack surface displacements averaged over a local volume 
(Allen, et al. 1987a). 

The presence of interply delamination in the laminate introduces jump 
discontinuities in the displacements and rotations of the normal line to the 
midplane. Since the normal to the midplane prior to deformation does not 
remain normal to the deformed surface of the midplane after deformation, it 
is no longer valid to neglect the out-of-plane shear strains in the problem for- 
mulation. In fact, it is through these out-of-plane shear strain terms that the 
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effects of the interply delaminations are introduced into the damage model 
development. The displacement equations are thus assumed to be of the fol- 
lowing form (Allen, el ah 1987c): 

M W) = u°{x.y) - :'{3° - 3(z - z x )tf) - H{z - c,) u ? , (2) 

v(z.y.z) = r°(i,y) - zrf - H(z - z^P) - H(z - z { )v? . (3) 

w(x.y.z) = u°(x.y ) - H(z - r,)u-, D . ( 4 ) 

where /, v ° . and w° are the undamaged midplane displacements: 3° and 
7 T are the undamaged ply rotations; up, i*P, and irp are the ply jump dis- 
placements due to dr. animation; 3® and 7 jP are the ply jump rotations due 
to del animations. The internai state variables corresponding to the delam- 
inauons are denned as the values of these jump rotations averaged over all 
the deiammation surfaces. Final]}*, H{z — z : ) is the Heavvside step function. 
The displacement equations are averaged over a local area to produce locallv 
averaged displacements that are used in the laminate formulation. 

The laminate constitution is then obtained by integrating the ply 
constitutive equations through the thickness. The ply constitution is assumed 
to be anisotropic since the jump discontinuities in the displacements resulting 
from delamination produce local anisotropic responses. That is, the local 
out-of-plane shear strains, and 7 i zzi resulting from delamination will 

contribute to the resultants. The laminate equations for the force and moment 
resultants are, 


{-'■} = MH4) - Wto} - {/"} + V D ), (5) 

{M} = [£]{<£} + \D]{k l ) + {»"> + {j D }, (6) 

where r 4i, [ B), and [ C) are the laminate extensional stiffness matrix, coupling 
stiffness matrix, and bending stiffness matrix, respectively. {/^} and {f®} 
are “damage induced forces’’ resulting from the matrix cracks and aelamina- 
tions. Their application to an undamage laminate will produce an equivalent 
amount of strain to those caused by the matrix crack and delamination dam- 
age surface kinematics, {g^} and {g D } corresponds to the “damage induced 
moments.’ 5 These “damage induced forces and moments” are determined from 
the corresponding internal state variables (Allen, et al. 1988). The laminate 
equations (5) and (6) are substituted into the plate equilibrium equations 
to yield the governing differential equations for the plate deformations. The 
governing equation for the out-of-plane deformation is (Buie 1988), 
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The in-plane governing equations are of similar form. The governing differ- 
ential equations are then integrated against a test function in the variational 
formulation of the laminated plate equilibrium equations. 

Finite Element Discretization 

In the current formulation, it is assumed that a total of five degrees of 
freedom exist at each node. The components of deformation at a node consist 
of two in-plane displacements, one out-of-plane displacement, and two out-of- 
plane rotations. The following displacement fields are assumed to represent 
the components of deformation within the element: 


”1 = !>>;'■ 

;=i 


( 8 ) 


vi = E (9) 

E*? 4 ?- < io) 

i=l 

where {£?} are the out-of-plane displacement and rotational components. ^ 
and represent the shape functions for the element, m is the number of nodes 
the element contains and p is three times the number of nodes. Substituting 
the aforementioned displacement fields into the weak formulation of the plate 
equilibrium equations will result in the following elemental stiffness and force 
matrices (Buie 1988), 




where [A - *] is the elemental stiffness matrix, {F 4 } is the applied force matrix, 
and {Am} and {F&} are the “damage induced” force matrices resulting from 
matrix cracking and delamination, respectively. For the model developed in 
the present work, a three node triangular element with five degrees of freedom 
per node is selected. This element is formed by combining the constant strain 
triangular element and a nonconforming plate bending element (Buie 1988). 

Internal State Variable Growth Law 

The values of the internal state variable used to represent the effects 
of matrix cracking and delamination are found through the damage evolution 
laws. These evolution laws describe the rate at which the internal state vari- 
ables are changing and are functions only of the current state at each material 
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point. The damage state, as described by the damage variables, can be found 
at any point in the loading history by integrating the damage evolutionary 
laws. The authors have recently developed one such relationship for the inter- 
nal state variable representing mode I matrix cracking under fatigue loading 
conditions (Lo, et al. 1990 y . The rate of cnange of this internal state variable 
during fatigue loading is expressed by 


do. 22 


da-2 

— kG^N 


( 12 ) 


where describes the change in the interna] state variable for a given change 
in the crack surface area, k and n are material parameters. G is the strain 
energy release rate for the damaged ply, and A is the number of load cycles. 
Because the evolutionary relationship is dependent only on ply level quantities, 
its application is not restricted to a particular laminate stacking sequence. The 
interactions with the adjacent plies and damage sites are implicitly reflected 
in the calculations through the laminate averaging process. Equation (12) can 
therefore be used to model both the transverse matrix cracks and axial splits 
in crossply laminates subjected to tensile loading conditions. The formulation 
of the delamination evolutionary relationship is currently in progress. It will 
be incorporated into the model once it becomes available. 

Model Application 

A finite element code has been developed based on the above formu- 
lation. The implementation of the damage evolutionary relationship shown in 
equation (12) requires that the solution algorithm be repeated for every load 
cycle. During each cycle, the ply stresses are calculated and used to determine 
the increment of change in the matrix crack internal state variable for each 
ply. The updated damage state is then used in the calculation of the laminate 
and ply response at the next load cycle. 

The response of a laminated tapered beam subjected to uniaxial fa- 
tigue loading conditions has been examined. This beam has a [ 0 / 9 0 2 ] ^ stack- 
ing sequence and possesses the material properties of continuouslv reinforced 
.454/3502 Graphite/Epoxv. Its length is 17.78cm and the width is 4.50cm at 
the clamped end and tapered to 2.87cm at the end where the load is applied. 
The beam is loaded at a distributed load amplitude of 17.5 x lO'N'/m and 
R = 0.1. Since the beam is symmetric about its length, it is sufficient to 
model half of the beam with a mesh containing 28 elements, 24 nodes, and 
120 degrees of freedom. The development of matrix cracks, during the loading 
history, in the 0° and 90° plies as represented by the internal state variables 
is shown in Figs. 1 and 2, respectively. For clarity, the damage levels in each 
ply are normalized by the largest value of the internal state variable within 
that ply at the end of 22495 load cycles. The predicted results indicate that 
the matrix cracks first occur in the 90° plies at the narrow end and progress 
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Fig. 1. The accumulation of matrix cracks in the 0° plies during various points 
in the loading history. 



Fig. 2. The accumulation of matrix cracks in the 90° plies during various points 
in the loading history. 
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toward the wider end as loading continued. Axial splits in the 0° plies do not 
develop until after the appearance of the transverse matrix crack in the 90° 
plies. Thus, the transverse matrix cracks serve as the driving mechanisms for 
the axial splits. As the transverse matrix cracks accumulate, loads carried by 
the 90° plies are transferred to the 0° plies. Depending on the amplitude of 
the latigue load, more than half of the load initially carried by the undamaged 
90° plies can b:* transferred to the 0° plies (Lo, et al. 1990). The additional 
loads coupled with the stress concentrations caused by the transverse matrix 
crac.Ks create suitable conditions for the growth of axial splits. The results 
also show that the progression of transverse matrix cracks alone the lenci;: 
of the tapered beam decelerates near the midway point. The stresses in the 
region beyond the mad way point are most likely to be insufficient in sustaining 
additional damage. Instead, the intensity of matrix damage increases at the 
narrow end during the latter portion of the loading historv. This process will 
continue until the matrix cracks have either reached the saturation level or 
when the laminate fails. The model predictions are qualitatively suDported 
by the experimental result shown in Fig. 3. 

Summary and Conclusion 

The authors have implemented a continuum damage mechanics model 
for laminated composite plates into a finite element computer code. The rep- 
resentation of matrix damage by second ordered tensor valued internal state 
variables and the utilization of modified laminate equations in the formulation 
has enabled the modeling of laminated composite plates with two dimensional 
rather than three dimensional finite elements. The reduction in the compu- 
tational requirements of this model over that of the three dimensional finite 
element analysis has made the analysis of composite plate structures subjected 
to complex loading histories a feasible task. The accumulation of matrix cracks 
in a tapered laminated composite beam loaded in fatigue is examined. The 
results reflect the damage induced load redistribution among the plies. Since 
this transfer of load influences the characteristic of damage accumulation and 
the eventual failure of the laminate, this model can assist in the determination 
of a laminated composite plate component's structural integrity. 
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